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Abstract 

We present the first fully quantum mechanical calculation of photon radiation 
accompanying charged particle decay from a barrier resonance. The soft- 
photon limit agrees with the classical results, but differences appear at next- 
to-leading-order. Under the conditions of o-decay of heavy nuclei, the main 
contribution to the photon emission stems from Coulomb acceleration and 
may be computed analytically. We find only a small contribution from the 
tunneling wave function under the barrier. 



1 



Typeset using REVI^X 



Nuclear fission and a-decay are interesting processes that involve both tunneling and 
the acceleration of charged particles in Coulomb fields. This raises the question of whether 
the tunneling process affects the bremsstrahlung emission. A semiclassical theory with an 
affirmative conclusion has been given by Dyakonov and Gornyi in ref . . In an experiment 
on the spontaneous fission of 252 Cf, Luke et. al. found a null result and gave an upper limit 
to the bremsstrahlung rate. In the case of a-decay of heavy nuclei, two recent experiments 
by D'Arrigo et. al. ||, and Kasagi et. al. [|J detected accompanying photon radiation. The 
latter authors claimed to observe interference effects with a tunneling contribution to the 
bremsstrahlung, interpreting their results in the framework of ref. . This gives urgency to 
carry out a full quantum mechanical calculation of the bremsstrahlung. We describe here 
a calculation with the following assumptions: a single-particle barrier model to describe 
the a-nucleus wave function; initial state treated as a Gamov state with a complex energy; 
perturbative coupling of the photon to the current, taken in the dipole approximation. 

It is convenient to use the acceleration form of the dipole operator, 

(fW) = ^r(f\[H,m = ^(f\^), (1) 

where H is the a-particle Hamiltonian with potential V and E 1 is the transition energy 
between initial and final state |z) and |/), respectively. The perturbative expression for 
dipole photon emission during decay may be obtained from Fermi's golden rule and is given 
by 

Here dP/dE^ is the branching ratio to decay with a photon emission, differential in the 
photon energy E 1 . The wave functions $«(r) and $/(?") are the radial wave functions of the 
initial and final state of the a-particle, respectively, with normalizations specified below. 
The effective charge Z e g for dipole transitions is given by Z c g = ((A — 4)z — A(Z — 2) /A 
where z = 2 is the charge of the a-particle and Z, A are the charge and mass number of the 
decaying nucleus, respectively. 



We take the potential in the single-particle Hamiltonian as the Coulomb outside a radius 
r and a constant inside, 

Zze 2 

V{r) = — e(r-r )-V e(r -r). (3) 

We will see later that the results are quite insensitive to the choice of parameters V and r , 

provided the decay properties are reproduced. The initial state $j is a resonant state of zero 

angular momentum normalized to a unit outgoing flux of particles. Its radial wave function 

i 

is given in terms of F and G Coulomb wave functions by 2 [Go(^, kr) + iFo(r] } kr)] jr 
outside r and is proportional to the jo( Kr ) spherical Bessel function inside. The Sommerfeld 
parameter 7] is given by rj = zZ r f™ it is much larger than one for heavy nuclei. The wave 



numbers k and k satisfy k = h~ 1 y/2mE a and k = h~ l ^2m(E a + Vo) where E a is the a- 
decay energy. Matching of the wave functions at r = r yields the amplitude of the inner 
wave function as well as the (complex) energy of the resonant state. 

The parameters r and V of our nuclear potential (|3|) are fixed to reproduce the empirical 
decay energy E a and mean life r of the decay. The mean life time depends on the parameters 
through the equation 

2E a r kr Gq(t}, kr ) f sin2/tr \ 



h 2 sin ktq 



kro 

In eq. @ and also in the wave function matching, we neglect terms with F Q and Fq which 
are of order 0(A) <^ 1 compared to G and G' . Here A = is a small parameter, and 
the primes denote derivatives with respect to kr. 

As is well known, there are multiple solution sets (r , Vo) for a given decay energy and 
mean life, distinguished by the number of nodes of the inner wave function ||. Typical 
solution sets for the nuclei of interest are shown in Table Q. Our simple model gives reasonable 
radii close to or slightly larger than the nuclear radius for Vo in the range to 150 MeV [Q . 
The results presented below do not depend on a specific choice of a solution. 

The continuum final states are normalized to give the unit operator when integrated 
over energy, 5(r — r') = J c/£'$£;(r)$£;(r / ). The radial wave function is a standing wave 



having the form ( J^P fc , j 2 (G\(rf , k'r) sin a + Fi(r]', k'r) cos a) jr outside r and proportional 
to the ji{n'r) spherical Bessel function inside. Primed quantities are defined similar to 
the corresponding unprimed quantities, but with the energy diminished by the emission 
of a photon. Matching the wave function at r = ro yields expressions for tana and the 
amplitude of the inner wave function. The final state is off resonance for almost all final 
energies (E a — E^), and thus tana ~ 0(A) <C 1, i.e. the final wave function does not 
penetrate the nucleus significantly and is a true continuum wave function. Outside the 
Coulomb-barrier the wave function is very well approximated by the regular Coulomb wave 
function only. 

The matrix element in eq. (0) can now be written down. It has a delta function contri- 
bution at r and an integral over the derivative of the Coulomb field outside, 

x [F^r]' , k'r ) + Gi(r]' , k'r ) tana] G (rj,kr ) 
-zZe 2 / drr- 2 { [Ft(rf, k'r) + Gi(V, k'r) tana] 

x [G ( V , kr) + iF ( V , kr)] } j + 0(A). (5) 

We separate the expression (|5|) into real and imaginary contributions and consider the latter 
first. 

We may neglect the contribution of the term F (i], kr)Gi(r]', k'r) tana to the integral 
since it is of order 0(A). Thus, the imaginary part is an integral over two Fj functions. 
Therefore it contains those contributions to the bremsstrahlung that stem from the classical 
acceleration in the Coulomb field. It can be treated analytically as follows. We first extend 
the lower limit of the integral to zero, which only introduces an error of the order 0(A). 
The resulting integral may be expressed in terms of hyper geometric functions as |],[7| 

POO Anf 

\ —F 1 {r ] ',k'r)F {T),kr) 
Jo r A 

= kk'{k'\l + ir)'\M - &|l + z77|Mi} (6) 

where 
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M = ( g y (r?+r?,) |r(j + i + V)||r(j + 1 + ^)| 

3 \V + V'J (k - k') 2 {2j + 1)! 

||| 2 F 1 fj + l-^,j + l-zV,2i + 2;-^V (7) 

Here 3F1 denotes the hypergeometric function and we have defined || 

£ = 7/ - 77. (8) 

In the limit of vanishing photon energy the imaginary part of the matrix element fl5|) may 
be computed directly, using [fTcR . 



lim Im(% l^yl = -J m f Q =2 - (9) 

The real part of the matrix element @ is a sum of two terms which, in contrast to the 
imaginary part, involve contributions from the irregular Coulomb wave functions Gj. Thus, 
it describes those contributions to the bremsstrahlung that are associated with tunneling. 



In the limit of vanishing photon energy, this amplitude reduces to [10.11 



fc Re <WI*W^7TTF (10) 

Notice that the dependence on the inner barrier parameters has disappeared. A comparison 
with the imaginary part (||) shows that the real part fllQl) is suppressed by a factor rj. For 
nonzero photon energy we have to treat the real part of the matrix element ([5]) numerically. 
However, the numerical evaluation shows that the real part still is suppressed in comparison 
to the imaginary part. This implies that only a smaller fraction of bremsstrahlung is emitted 
during tunneling. Note also that the contributions associated with classical acceleration and 
tunneling do not interfere since they differ in phase by i. 

We will now make the connection to semiclassical and classical limits. For heavy nu- 
clei, the Sommerfeld parameters r] are large and the Coulomb wave functions Fj may be 
approximated by their WKB-wave functions 

F^ B ( V ,kr) = (k 2 /f(r)f sin <P, (11) 
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with 



f(r) = k 2 -2kr]/r - j(j + l)/r 2 and (12) 



7T 



4 +l /t dr>lHr>)]- (13) 



In leading order in rj, rj' one finds 

dr 



J^Fr B (v\k'r)Fr B (v,kr) 



-e 2 5 



k + k' rj 



(14) 



where e = - + ^ 4 ^ and 77 = (rf + rj)/2. K v {z) denotes the modified Bessel function and K' v 
its derivative with respect to the argument. 

A comparison of the semiclassically evaluated integral fll4l) with the quantum mechani- 
cal result (|6]) shows that they deviate from each other by less than one percent for photon 
energies E 7 up to 1 MeV. We recall that the semiclassical computation neglects any contri- 
butions from the wave functions at radii smaller than the classical turning point, i.e. any 
contribution from the tunneling. This clearly justifies the attribution of tunneling to the 
real part of the matrix element, (^j), alone. 

Next we consider the classical and the soft photon limit. The classical formula valid at 



all frequencies can be derived from [12 



dP 2aZ 2 s \I(cu)\ 2 



dE^ 3n E^y 
with / the Fourier transform of the time-dependent acceleration, 

f°° , dv 



(15) 



POO (£y 

I(uj) = c~ 1 dt — exp(icut). (16) 
Jo dt 



This integral can be expressed in terms of the dimensionless parameter 



as 



= ] J^I dzexp ^ (r=v + artanh z ) • ( 18 ) 



In the limit of small photon energy we find 

dt^ = AaZl L E^ 

dE 1 3tt mc 2 7 ' 1 ' 

Because this only depends on the asymptotic motion of the particles, the quantum result 
must coincide. Inserting the results @, fllOD in eq. @ indeed yields the classical result, eq. 



(ED- 

More interesting is to examine the next-to-leading ^-dependence and compare the 
quantum and classical behavior. It turns out that the classical parameter ( in eq. ([[7]) is 
essentially the same as the quantum small parameter £ defined in eq. (§) (( = 2£ + 0(£ 2 )). 
This parameter may also be identified with the product of the photon frequency and the 
barrier tunneling time. The (^-dependence of the classical and quantum calculations are 
compared in Fig. [1|. The solid line shows the classical prediction (|T5|). The dashed and 
dotted lines show the quantum result with and without tunneling contributions, respectively. 
We see that the tunneling contributions remains small even at a finite photon energy. 

One might have expected that the classical curve would be tangent to the quantum at 
uj = 0: in scattering bremsstrahlung is determined by on-shell amplitudes to next-to-leading 
order fl3fl . We find that the two curves are indeed very close in the neighborhood uj = 0, 



but the slopes are not identical. For large photon energies the classical result overestimates 
the photon emission rate considerably since the classical formula (pl|) neglects any energy 
loss of the escaping a-particle. This point has been discussed in the framework of photon 
emission in spontaneous fission by Luke et. al. 0, and earlier in the framework of Coulomb 
excitation by Alder et. al. [0. 

The quantum mechanical results for 214 Po and 226 Ra are practically identical to those for 
210 Po when plotted as in Fig. [1], normalized to the uj = rate (|1^) and plotted as a function 
of (. Since ( is inversely proportional to the decay energy E a , the rates are higher for higher 
decay energies. Thus for 214 Po decay, with an a-decay energy of 7.7 MeV, the predicted 



rate for E 1 = 0.6 MeV is 65 times higher than for 210 Po. 

Finally, we compare the results obtained in this work with experiment. In the case of 
210 Po, our result displayed in Figure is consistent with the experimental result obtained 
by Kasagi et. al. [|J suggesting that no interference resulting from photon emission during 
tunneling is needed for an explanation of the experiment. In the case of 214 Po and 226 Ra, 
D'Arrigo et. al. [p| reported photon emission rates that are larger than expected from the 
classical formula (|i~5f) . Thus, their results are also more than one order of magnitude larger 
than our quantum mechanically computed rate. We cannot trace the origin of this difference. 

In summary, we have used Fermi's golden rule to compute the emission of bremsstrahlung 
in a-decay of heavy nuclei. The dominant contribution to the photon emission rate 
stems from classical acceleration and is given in closed form. Only a smaller fraction of 
bremsstrahlung is emitted during tunneling. This finding is consistent with experimental 
data on 210 Po. 

We thank R. Vandenbosch, A. Bulgac and N. Takigawa for discussions. This work was 
supported by the Department of Energy on contract No. DE-FG-06-90ER-40561. 
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TABLES 





214p 


210p o 


226 Ra 


E/MeV 


7.7 


5.3 


4.8 


A = h/rE a 


3.7- i(r 19 


7.3 • i(r 30 


1.9 • 10~ 33 


ro/fm 


9.19 


8.76 


9.75 


V /MeV 


12.1 


16.7 


12.9 



TABLE I. Parameters of the nuclear potential. The table lists the parameters ro and Vq of 
the nuclear potential (||) that were used in the present calculations. For each listed nuclei, the 
presented parameters reproduce the experimentally known values of the a-particle energy E a and 
the mean life time r. 
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FIGURES 
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FIG. 1. Comparison of classical and quantum mechanical photon emission probability in 
a-decay of 210 Po. Curves show the probabilities normalized to the low-energy formula, eq. (|l9|). 
as a function of the scaled photon frequency £ defined in eq. (|l~7|). The classical probability, eq. 
(|l5|), is shown as the solid line. The quantum probabilities (nearly exponentially falling lines) are 
shown for the full quantum mechanical treatment (dashed) and for the approximation that neglects 
contribution from tunneling (dotted). £ = 1 corresponds to E~ ~ 0.24MeV. 
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FIG. 2. Photon emission probability comparing the quantum calculation with experiment. 
Logarithmic plot of -gg- as a function of photon energy for 210 Po (full line). The experimental 
data for 210 Po (datapoints with errorbars) are taken from Ref. Jl. 
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